Surface Reconstructions and Bonding via the Electron Localization Function: 

The Case of Si(OOl) 
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The bonding pattern of a covalent semiconductor is disrupted when a surface is cut while keeping 
a rigid (truncated bulk) geometry. The covalent bonds are partly reformed (with a sizeable energy 
gain) when reconstruction is allowed. We show that the "electron localization function" (ELF) — 
applied within a first-principles pseudopotential framework — provides un unprecedented insight into 
the bonding mechanisms. In the unreconstructed surface one detects a partly metallic character, 
which disappears upon reconstruction. In the surface reformed bonds, the ELF sharply visualizes 
strongly paired electrons, similar in character to those of the bulk bonds. 
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I. INTRODUCTION 

The driving mechanism for surface reconstructions of a 
covalent material is the tendency of low-coordinated sur- 
face atoms to saturate dangling bonds. It is then obvious 
that a completely different bonding pattern differenti- 
ates the unreconstructed (i.e. bulk terminated) surface 
from the reconstructed one. When several different re- 
constructions are compared, these show in turn different 
bonding features. First-principles electronic structure 
calculations have provided over the years a major insight 
into these phenomena: however, the major tools used in 
investigating bonding at large has been the analysis of 
either the electronic charge density or the projected den- 
sity of states. Since a few years a much more informative 
tool has been introduced in the quantum chemical litera- 
ture to deal in a quantitative way with bonding features: 
this goes under the name of "electron localization func- 
tion" (ELF). 1-3 Although bonding at a reconstructed 
surface looks like an ideal arena for an ELF investiga- 
tion, we are aware of only one such study, which how- 
ever is performed at the non-selfconsistent tight -binding 
level. 4 We present here a state-of-the-art thorough inves- 
tigation about a paradigmatic case: namely, the silicon 
(001) surface. We show that the insight into the na- 
ture of bonding at the surface is very accurate and clear 
when the ELF is analyzed. We demonstrate — through 
a series of contour plots in high-symmetry planes — the 
outstanding ability of ELF in discriminating between 
the surface bonds occurring upon dimerizations at the 
Si(001) surface, thus providing an unprecedented insight 
into the physical mechanisms which drive the reconstruc- 
tion. Dangling orbitals and surface bonds are visualized 
with a resolution incomparably sharper than by using the 
current tools, such as charge-density plots or projected 
densities of states. 

We are using here the approach which has become 
the "standard model" in first-principle studies of cova- 
lent semiconductors: 5 namely, density-functional theory 
(DFT) 6 ~ 8 with norm-conserving pseudopotentials. 9 This 



approach has got a dominant role in the literature be- 
cause of many of its key features, but ELF investigations 
within it have been quite scarce so far: ELF has been 
originally proposed as an all-electron tool, and as such 
has been mostly applied. Instead, we are going to show 
in the present work that by getting rid of the core elec- 
trons and focussing on the bonding electrons only, one 
makes the ELF message particularly perspicuous. 

We study here an elemental system, where the key is- 
sue is metallic versus covalent bonding. Since the early 
days of electronic-structure theory, we understand bond- 
ing in simple metals through the paradigm of pseudopo- 
tential perturbation theory: 10 ' 11 the valence electrons be- 
have basically as a free-electron gas in the region out- 
side the ionic cores. At the opposite extreme, we under- 
stand covalent bonding through the paradigm of the H2 
molecule, whose electrons are strongly paired in a singlet 
state; in a more general case, covalent bond is character- 
ized by "localization" of valence-electron pairs in appro- 
priate regions of space. While charge-density plots do 
not help much in discriminating metallic from covalent, 
the ELF provides a quantitative measure of "metallic- 
ity" vs. "covalency" of a given bond, or more generally 
of a given valence region of the system. This measure is 
particularly significant in a pseudopotential framework 
where — amongst other things — a simple metal immedi- 
ately displays the free-electron nature of its valence elec- 
trons. We furthermore observe that ELF is by definition 
an orbital-free approach: this feature is preserved when a 
first-principle framework (as opposed to a tight -binding 
one) is adopted, as it is done here. 

The case study chosen for this work is possibly the the- 
oretically best known semiconductor surface. Its most 
relevant features are therefore very well understood in 
the literature: in particular, the unreconstructed surface 
is metallic and the reconstructed one is insulating. 12 The 
aim of the present paper is not to demonstrate anything 
at variance with the common wisdom. Instead, our aim 
is to show how the known features of this surface can be 
recovered in a simple and meaningful way from an ELF 
analysis. It is also worth to stress that the ELF does 
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not make use of the spectrum of the system, while in- 
stead is a pure ground state property. In perspective, we 
are proposing ELF, in a modern pseudopotential frame- 
work, 5,9 as a powerful tool for investigating bonding be- 
tween over- or under-coordinated atoms in more com- 
plex, and worse known, situations. 



II. ELECTRON LOCALIZATION FUNCTION 

ELF was originally introduced by Becke and Edge- 
combe as a measure of the conditional probability of find- 
ing an electron in the neighborhood of another electron 
with the same spin. Starting from the short-range be- 
havior of the parallel-spin pair probability, they defined 
a scalar function, conveniently ranging from zero to one, 
that uniquely identifies regions of space where the elec- 
trons are well localized, as occurs in bonding pairs or 
lone-electron pairs. By definition, ELF is identically one 
either in any single electron wavefunction or in any two- 
electron singlet wavefunction: in both cases the Pauli 
principle is ineffective, and the ground wavefunction is 
nodeless or, loosely speaking, "bosonic" . In a many- 
electron system ELF is close to one in the regions where 
electrons are paired to form a covalent bond, while is 
small in low-density regions; ELF is also close to one 
where the unpaired lone electron of a dangling bond is 
localised. Furthermore, since in the homogeneous elec- 
tron gas ELF equals 0.5 at any density, values of this 
order in inhomogencous systems indicate regions where 
the bonding has a metallic character. 

A major advance is the ELF reformulation due to Savin 
et al., 13 which is adopted in the present work. Owing to 
the Pauli principle, the ground-state kinetic energy den- 
sity of a system of fermions is no smaller than the one 
of a system of bosons at the same density: ELF can be 
equivalently expressed in terms of the Pauli excess en- 
ergy density, thus making no explicit reference to the 
pair distribution. This is particularly convenient in our 
case, since the pair density is outside the scope of DFT. 
The Savin et al. reformulation also provides a meaningful 
physical interpretation: where ELF is close to its upper 
bound, electrons are strongly paired and the electron dis- 
tribution has a local "bosonic" character. 

The definition of ELF is: 



£(v) = 



1 



l + [D(v)/D h (v) 



(1) 



thus taking by design values between zero and one. Fol- 
lowing Savin et al., D(r) is the Pauli excess kinetic en- 
ergy density, defined as the difference between the kinetic 
energy density and the so-called von Weizsacker kinetic 
energy functional: 14 



D(r) = -V r V, 



'P(r, r') 



Vn(r)| : 
n(r) 



(2) 



where p is the one-body reduced (spin-integrated) den- 
sity matrix. The von Weizsacker functional provides 



a rigorous lower bound for the exact kinetic energy 
density 1 and is ordinarily indicated as the "bosonic" 
kinetic energy, since it coincides with the ground-state 
kinetic energy density of a non-interacting system of 
bosons at density n(r). Therefore, D is positive semidef- 
inite and provides a direct measure of the local effect of 
the Pauli principle. The other ingredient of Eq. (1) is 
Dh(r), defined as the kinetic energy density of the ho- 
mogeneous electron gas at a density equal to the local 
density: 



D h (r) = A(37r 2 )ln(r)l. 



(3) 



It is obvious to verify that the function £ is identical to 
one in the ground state of any two-electron system, while 
is identical to 0.5 in the homogeneous electron gas at any 
density. 

As commonly done in many circumstances — including 
in other ELF investigations 13 — we approximate the ki- 
netic energy of the interacting electron system with the 
one of the nonintcracting Kohn Sham (KS) one. We 
therefore use the KS density matrix: 



p(r,r') = 2^^(r)^ s (r'), 



(4) 



where <j>f s (r) are the occupied KS orbitals. Such ap- 
proximation is expected to become significantly inaccu- 
rate only in the case of highly correlated materials. 

The original ELF definition is an all-electron one, and 
has the remarkable feature of naturally revealing the en- 
tire shell structure for heavy atoms. Such a feature is of 
no interest here, since only one electronic shell (the sp 
valence one) is involved in the bonding of Si atoms in 
any circumstances. The pseudopotential scheme simpli- 
fies the landscape, since only the electrons of the relevant 
valence shell are dealt with explicitly: the ELF message 
comes out therefore much clearer, with basically no loss 
of information. 15 In the spatial regions occupied by core 
electrons, the pseudo-electronic distribution shows a de- 
pletion and ELF assumes very low values. Ouside the 
ionic cores, in the regions relevant to chemical bonding, 
the norm conservation endows the pseudocharge density 
with physical meaning, as widely discussed in the mod- 
ern pseudopotential literature: 5 ' 9 in the specific case of Si, 
there is an experience of about 20 years in dealing with a 
plethora of different physical problems. By the same to- 
ken of norm conservation, even the pseudo-orbitals — and 
hence their kinetic energy density — closely map the all- 
electron ones in the bonding regions. The pseudo ELF 
carries therefore — in the material of interest here — the 
same information as the all-electron ELF, while it re- 
moves irrelevant and confusing features due to the inner, 
chemically inert, shells. 
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III. CALCULATIONS 

The (001) silicon surface is strongly reconstructed, 
with several different reconstructed structures. A rela- 
tively simple (2x1) reconstruction is formed by surface 
atoms moving together in pairs to form dimers. 




Q 1st layer O 2nd layer 

• 3rd layer • 4th layer 

FIG. 1. Side (left) and top (right) views of the (2 x 1) recon- 
structions of the Si(001) surface, (a) Buckled-dimer, (b) sym- 
metric-dimer, and (c) ideal (bulk-terminated) surface. The 
atoms in the first four layers are shown, with symbols as in- 
dicated in the figure bottom. The surface unit cells (shaded 
areas) are also shown. 

The driving force for the dimer formation is the elimi- 
nation of a dangling bond from the surface atoms. Each 
atom of the unreconstructed surface is bonded to only 
two neighbours and therefore has two dangling bonds 
projecting out of the surface. If the atoms move together 
in pairs, forming a new bond between them, then one of 
these dangling bonds will be eliminated from each mem- 
ber of the pair. This leads to a considerable energy gain 
since half of the broken surface bonds are reconstructed. 
The buckled dimer model is a simple modification of the 
dimer model in which each dimer is tilted or "buckled" 
out of the plane of the surface. We have considered there- 
fore the unreconstructed and the (2 x 1) symmetric and 
buckled dimer surfaces, in order to better understand the 
connection between surface reconstruction and electronic 
localization. A schematic side and top view of these three 
surfaces is shown in Fig. 1. 
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FIG. 2. ELF contour plot for the (001) bulk terminated 
silicon surface, in the plane containing two nearest-neighbor 
atoms: the relative dangling orbitals are revealed by the 
black ELF regions. The grey-scale is also shown: dark 
(clear) regions correspond to large (small) ELF values. In the 
schematic side view, black balls and solid lines correspond, 
respectively, to in-plane atoms and bonds; white balls and 
dashed lines to off-plane atoms and bonds. 

All the calculation use quite standard ingredients: a 
plane-wave expansion of the KS orbitals with a 10 Ry 
kinetic-energy cut-off, 16 k-points on a Monkhorst and 
Pack 16 mesh for the irreducible Brillouin zone integra- 
tion, and a norm-conserving pseudopotential in fully 
non-local form. 17 All our calculations were performed in 
a supercell geometry where the surface is modeled by a 
finite-size slab periodically repeated in the direction nor- 
mal to the surface. For the ideal (bulk terminated) sur- 
face we use a unit cell containing 11 layers of atoms sep- 
arated by the equivalent of 13 atomic layers of vacuum. 
For the symmetric and buckled dimer (2 x 1) surface, we 
have reproduced the Roberts and Needs 18 calculation, 
using directly their relaxed atomic coordinates. 




FIG. 3. Contour plots for bulk aluminum in the [100] 
cristalline plane. Left panel: ELF, where the same grey-scale 
as in all ELF figures is adopted. Right panel: pseudocharge 
density. 
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FIG. 4. Pseudocharge density contour plot for bulk silicon 
in the same crystallographic plane as for Fig. 2. 

We show in Fig. 2 the ELF contour plot for the unre- 
constructed surface, in the plane containing the dangling 
orbitals and two nearest neighbor surface atoms. As al- 
ready discussed above, pseudoatoms correspond to the 
white regions (ELF minima). Before discussing the sur- 
face results, we pause to illustrate the bulk features and 
the capability of ELF in discriminating between metallic 
and covalent bonding. 




FIG. 5. ELF contour plot for the (001) bulk terminated sil- 
icon surface. Plane orthogonal to the dangling bond, passing 
through the center of the surface unit cell. For the schematic 
side view, symbols as in Fig. 2. 

In the bulk Si region, one clearly sees the almost 
black regions in Fig. 2 between nearest-neighbor atoms, 
whose bonding pattern in this geometry has the shape 
of a "zig-zag" chain. In these bonding regions our cal- 
culated ELF attains the maximum value of 0.96, thus 
indicating that the Pauli principle has little effect. In 
agreement with a chemical picture of the covalent bond, 
we associate these regions to the opposite-spin electron 
pairs — actually a "bosonic" system — localized between 
every pair of bonded atoms. This is to be contrasted to 
metallic bonding, where the valence electrons have a free- 
electron nature: to make this point better clear, we take 
as an example a paradigmatic simple metal: crystalline 
aluminum, whose bulk ELF is shown in the left panel 
of Fig. 3 to the same grey scale. Comparing this plot 



with the bulk region of silicon in Fig. 2 we immediately 
appreciate the spectacular ELF ability to distinguish in 
a very clearcut way between metallic bonding and cova- 
lent bonding. The ELF plot in aluminum shows — outside 
the core regions — a large grey area, which correspond to 
a jcllium-like (or Thomas-Fermi) ELF value. Actually, 
the maximum value attained by £(r) between nearest — 
neighbor atoms is only 0.61. 

Comparison of the two ELF plots provides therefore 
the most significant and perspicuous visualization of 
the important qualitative difference between the covalent 
bond and the metallic one. In the bulk of the two mate- 
rials the Pauli principle plays quite a different role. The 
other typical tools for analysis, such as charge-density 
plots or projected density of states, lack by far a similar 
sharpness. As an example, we show the corresponding 
pseudocharge density plots in the right panel of Fig. 3 
(aluminum) and in Fig. 4 (silicon). From such figures, 
hardly any information about metallicity or covalcncy 
can be drawn. 



IV. RESULTS 

We now move to discuss the surface results. Focussing 
on the surface region in Fig. 2, one notices that the two 
dangling orbitals are also well visualized by the ELF plot: 
every dangling bond is associated to an isolated elec- 
tron, thus corresponding to an high ELF value. How- 
ever, a more striking bonding feature of the same unre- 
constructed Si(001) surface is detected when analyzing 
ELF in the orthogonal plane, passing midway between 
two nearest neighbor surface atoms (Fig. 5). There is 
a grey "fermionic channel" — defined as a region where 
ELF is of the order of 0.5 — between the two nonbonded 
atoms: this channel in the midpoint has a perfect circu- 
lar section, but it is actually formed by two grey-ELF 
strips going round the surface atoms and intersecting at 
the midpoint, as clearly shown by a top view (left panel 
of Fig. 6). 




FIG. 6. ELF contour plot for the Si(001) bulk terminated 
(left panel) and symmetric-dimer reconstructed (right panel) . 
Surface plane passing through the topmost atoms. 

The presence of this extended metallic-like system 
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along the unrelaxed surface must be connected to the 
strong surface strain, due to the nonbonded atoms. As 
soon as the surface is relaxed and the dimer formation is 
allowed, there is a strong surface "bosonization" , associ- 
ated with a sizeable reduction of the fermionic channel. 
This localization effect due to the dimerization process is 
perspicuously shown in Fig. 6, where the unreconstructed 
situation is compared to the reconstructed one, for the 
symmetric-dimer case. The metallic character of the 
unreconstructed surface, and the insulating character of 
the reconstructed one, are well documented in the liter- 
ature. 12 What is remarkable here is that ELF visualizes 
such characters without using any spectral information 
about the system. 
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FIG. 7. ELF contour plot for the Si(OOl) (2 x 1) symmet- 
ric-dimer surface. Plane orthogonal to the surface, containing 
the dimer. For the schematic side view, symbols as in Fig. 2. 

A more detailed analysis of the ELF variations induced 
by rebonding at the surface is provided by Fig. 7 for 
the symmetric-dimer case and in Fig. 8 for the buckled- 
dimer one, both drawn in the plane passing through the 
surface atoms, and to be compared to the analogous plot 
of Fig. 2 for the unreconstructed case. In both cases the 
dimer has a strongly covalent bond character and is sur- 
rounded by a region of high ELF values: such localization 
effect is clearly connected to the surface reduced coordi- 
nation, since the ELF around the isolated topmost atoms 
has a more pronounced atomic character. 19 





FIG. 8. ELF contour plot for the Si(OOl) (2 x 1) buckled 
dimer surface. Plane orthogonal to the surface, containing 
the dimer. For the schematic side view, symbols as in Fig. 2. 

In conclusion, we have shown how the electronic sys- 
tem at a paradigmatic semiconductor surface changes 
its character upon reconstruction. ELF clearly visual- 
izes how electron pairing accompanies rebonding. In the 
unreconstructed surface, besides the dangling-bond elec- 
trons, we have detected channels of metallic-like charac- 
ter in the surface plane. Both features disappear upon 
reconstruction. 
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